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Abstract. Bianchi VIII vacuum solutions to Einstein's equations are causally 
geodesically complete to the future, given an appropriate time orientation, and 
the objective of this article is to analyze the asymptotic behaviour of solutions 
in this time direction. For the Bianchi class A spacetimes, there is a formu- 
lation of the field equations that was presented in an article by Wainwright 
and Hsu, and we will analyze the asymptotic behaviour of solutions in these 
variables. We also try to give the analytic results a geometric interpretation by 
analyzing how a normalized version of the Riemannian metric on the spatial 
hypersurfaces of homogeneity evolves. 



1. Introduction 

Let us begin by defining what we mean by Bianchi class A spacetimes. We do 
this by defining the concept class A vacuum initial data. Let us introduce some 
terminology. Let G be a 3-dimensional Lie group, e^, i = 1,2,3 be a basis of the 
Lie algebra with structure constants determined by [e^, ej] — lijSk- If lik — 0; then 
the Lie algebra and Lie group are said to be of class A and 

(1) 7^- = e.,™^'™ 
where the symmetric matrix n*^ is given by 

(2) „».^l^(;e^)«. 

Definition 1. Class A vacuum initial data for Einstein's equations consist of the 
following. A Lie group G of class A, a left invariant metric g on G and a left 
invariant symmetric covariant two-tensor fc on G satisfying 

(3) Rg-kijk'^ + (tTgk)^ =0 

and 

(4) V^trgfc - V^fcy = 

where V is the Levi-Civita connection of g and Rg is the corresponding scalar 
curvature, indices are raised and lowered by g. 

Consider class A vacuum initial data. We can choose a left invariant orthonormal 
basis {ci} with respect to g so that the corresponding matrix n*-' defined in (||) 
is diagonal with diagonal elements ni, n2 and n^. By an appropriate choice of 
orthonormal basis ni,n2,n3 can be assumed to belong to one and only one of the 
types given in table ^ We assign a Bianchi type to the initial data accordingly. 
This division constitutes a classification of the class A Lie algebras, see Lemma ^. 
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Table 1 



Bianchi class A. 



Type 
I 



ni 

+ 





n2 



+ 
+ 
+ 
+ 








II 



VIo 
VIIo 
VIII 
IX 



+ 



+ 
+ 
+ 



We can divide the initial data into three classes, 

1. Bianchi IX initial data. In this case the maximal globally hyperbolic develop- 
ment of the initial data is future and past causally geodesically incomplete. It 
expands for a certain period of time, reaches a moment of maximal expansion 
and then recoUapses. 

2. Bianchi I and VIIq initial data with tr^fc = 0. The corresponding develop- 
ments are causally geodesically complete. 

3. The remaining types of initial data. With a suitable time orientation such 
data evolve into a future causally geodesically complete spacetime expanding 
forever into the future. The development is however past causally geodesically 
incomplete. 

In this paper we are interested in the third class of initial data. 

There is a formulation of Einstein's equations due to Ellis and MacCallum [Q cov- 
ering, among other things, the Bianchi class A spacetimes. In an article by Wain- 
wright and Hsu a normalized version of the corresponding variables together 
with a different time coordinate were introduced. The main part of this paper 
consists of an analysis of the asymptotic behaviour of these variables. Our main in- 
terest is in Bianchi VIII, but we also consider the other Bianchi class A spacetimes 
in the third class as above. The Bianchi I, II and VIq cases are quickly handled, and 
are included for completeness, but it takes some time to analyze the asymptotics 
of Bianchi VIIq. Observe that Bianchi VIIq perfect fluids have been considered in 
H). Bianchi VIII is dealt with definitively, in the sense that it is proven to what 
the variables of Wainwright and Hsu converge. However, in certain situations this 
may not necessarily be enough. For instance, it is proven that N2 + N3 converge 
to infinity, if A^2, N3 > and iVi < 0, and that N2 — N3 converge to zero. It is 
conceivable that if one is interested in computing curvature invariants for instance, 
one has to consider expressions of the form N2 — , concerning which the results 
contained in this paper say nothing. 

Let us also mention that even though it does not require much time to sort out 
the asymptotics for Bianchi VIq, the results yield an interesting consequence. In 
fact, one gets the conclusion that Bianchi VIII is very different from Bianchi IX as 
far as the behaviour towards the singularity is concerned. In Q it was proven that 
generic Bianchi IX solutions converge to an attractor, where by the attractor, we 
mean the set of vacuum type I and II points, and the generic points are obtained by 
subtracting a finite number of positive codimension submanifolds from the manifold 
of Bianchi type IX points. It was also proven that for generic Bianchi IX solutions, 
the convergence to the attractor is almost monotone in the following sense. Given 
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an e > 0, there is a (5 > such that if x constitutes Bianchi IX initial data closer 
to the attractor than 5, then applying the flow to x in the past time direction 
will not result in points further away from the attractor than e. By distance, we 
here mean the ordinary Euclidean distance in the Wainwright Hsu variables in M^, 
the relevant space in the Bianchi IX matter case. This result, together with the 
result that generic Bianchi IX solutions have an a-limit point on the Kasner circle, 
yield the conclusion that generic solutions converge to the attractor. The almost 
monotone convergence constitutes the most difficult step. In Proposition |^, we 
prove that the corresponding statement for Bianchi VIII is not true. In fact, we 
prove that there is a fixed number 77 > 0, such that for any e > we can find 
Bianchi VIII initial data closer to the Kasner circle than e which will at one point 
to the past be at a distance further than rj away from the attractor. This does of 
course not prove that Bianchi VIII solutions may not converge to the attractor, but 
it does imply that the main part of the argument concerning Bianchi IX presented 
in is of no help in proving convergence to the attractor, and that a completely 
different method is required, assuming one still believes that generic Bianchi VIII 
solutions converge to the attractor towards the singularity. 

Let us try to explain the argument concerning the future asymptotics of a Bianchi 
VIII solution with iVi < and N3 > 0. By a monotonicity argument, one 
can prove that A^i — )■ and -^3 ^ 00 as r ^ 00. This is a very important 
observation in that it introduces a natural concept of order of magnitude; powers 
of N2 + N3. Writing the constraint as 

S+ + + ^[Nf + {N2 - m f - 2Ni{N2 + N3)] = 1, 

we conclude that all the terms involved are bounded since they are all non-negative. 
The derivatives of S+ and Ni{N2 + N3) are bounded by constants independent of 
the initial data, but we have 

i' = -3(iV2 + N3)y + y' = 3(7V2 + iVg)^ + ... 

where x — S_ and y — \/3{N2 — A^3)/2, and the dots represent expressions that 
can be bounded by constants independent of which Bianchi VIII solution we are 
considering. Thus, if (x^ + y'^Y^^ is big in comparison with {N2 + A^a)"^, y and x 
will essentially behave as sine and cosine, and the frequency of the oscillation will 
go to infinity as r ^ 00. The oscillations in x and y are not so interesting, but 
the evolution in S+ and Ni{N2 + N3) is. For that reason, it is natural to consider 
what happens from, say, one time at which I]_ = to the next. Before doing so, 
some effort has to be put into estimating the error committed in approximating y 
and X with sine and cosine, into proving that there actually are zeros of E- and 
estimating the variation of different objects during the time that elapses between 
the zeros. After having done that, it will be possible to express S+ and Ni{N2 + N3) 
at a zero of E_ in terms of the values of the same expressions in the previous zero. 
In other words, it will be possible to replace the continuous flow with a discrete 
map depending on only two variables. This approximation will not always be valid, 
but when it is not, it can be replaced by less sophisticated arguments. The discrete 
map is then the main tool in a sequence of technical arguments proving the desired 
result. 

When we evolve initial data we get a Lorentz metric g — —dt^ +9{t) on / x G where 
/ is an open interval and g(t) is a Riemannian left invariant metric on G. Let F be a 
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subgroup of the group of diffeomorphisms of G, acting properly discontinuously on 
G. We can consider F as acting on / x G by letting T act only on G. In situations 
where g(t) is invariant under F, we can take the quotient to get a solution to 
Einstein's equations on / x M, where M = G/F. Assuming M to be compact we 
can define the reduced Hamiltonian, c.f. 

iJroduccd = -{trgk)\o\{M,g), 

where k{t) is the second fundamental form of {t} x M. We can consider it to be 
a function of t. This object is of interest in the context of the work by Fischer 
and Moncrief ||^. We prove that the reduced Hamiltonian converges to zero for 
all the initial data in class three, though it should be remarked that the spatial 
hypersurfaces should be of Yamabe type -1, see for this result to be of interest. 
Finally, let us refer the reader to [|| and references therein for further discussion of 
Bianchi VIII vacuum solutions. 

The equations of Wainwright and Hsu are to be found in section ^ We mention 
some properties and describe an important tool in the analysis of the asymptotics of 
solutions to these equations; the monotonicity principle. In section ^ we relate the 
Lorentz metric on the development with the Wainwright Hsu variables. A formula 
expressing the reduced Hamiltonian in terms of these variables is also given. The 
main part of this paper consists of an analysis of the asymptotic behaviour of 
solutions to the equations of Wainwright and Hsu. This analysis occupies the 
remaining sections except for the last. The final section contains an analysis of 
the asymptotic behaviour of the Ricmannian metric on the spatial hypersurfaces of 
homogeneity. 



2. The equations of Wainwright and Hsu 



We formulate the equations we will use here and state some properties. The follow- 
ing section contains a derivation. Einstein's vacuum equations take the following 
form in the formulation due to Wainwright and Hsu 

N[ = (g-4E+)iVi 

= (g + 2S+ + 2\/3I]_)A^2 

(5) A^^ = (g + 2E+-2\/3I]_)A^3 
K = -{2-q)^+-3S+ 

T.'_ = -(2- - SS*- 
where the prime denotes derivative with respect to r and 
q = 2(S^ + E^) 

(6) S+ = ^[{N2 - Nsf ~ Ni{2N, ~ N2 - Ns)] 

S- = ^(A^3- A^2)(A^i-iV2- A3). 
The vacuum Hamiltonian constraint is 

(7) + + ^[a2 + A| + - 2(A^iA^2 + N2N3 + N1N3)] = 1. 

The above equations have certain symmetries described in Wainwright and Hsu Q . 
By permuting Ai, A^2, A^3 arbitrarily we get new solutions if we at the same time 
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carry out appropriate combinations of rotations by integer multiples of 27r/3 and 
reflections in the S-)-plane. Changing the sign of all the Ni at the same time 
does not change the equations. Classify points (A^i, -/V2, -^3, S+, £_) according to 
the values of TVi, A^2, -^3 in the same way as in table |l[ Since the sets Ni > 0, Ni < 
and Ni = are invariant under the flow of the equation we may classify solutions 
to (||)-(0) accordingly. When we speak of Bianchi VIII solutions we will assume 
two Ni > and one < 0. The Kasner circle is the subset of E_|_E_iViiV2iV3-space 
where N = and T,l + Y.'i = 1. 

We only consider solutions to (||)-(0) which are not of Bianchi IX type. By the 
constraint (Q) we conclude that q <2 for the entire solution. As a consequence the 
Ni cannot grow faster than exponentially due to (||) . The solutions we consider can 
thus not blow up in a finite time so that we have existence intervals of the form 
(— c», 00). 

Definition 2. Let /Ci, IC2 and /C3 be defined as the subset of the Kasner circle on 
which q - 4E+ <0, q + 2S+ + 2 vll]_ < and q + 2S+ - 2^31;- < respectively, 
c.f. (D. 

The set Y,^ = 0, N2 — N^ is invariant under the flow of Applying the 

symmetries to this set we get new invariant sets. 

The concepts a- and w-limit set will be useful. 

Definition 3. Let / G C°°(R",IR"), and consider a solution x to the equation 

dx 

(8) — ^ f ox, x{0) = xo, 

with maximal existence interval We call a point x^ an Lu-limit point of the 

solution x, if there is a sequence tk t+ with x{tk) — > x*. The uj-limit set of x 
is the set of its w-limit points. The a-limit set is defined similarly by replacing t+ 
with 

Lemma 1. Let f and x he as in Definition^. The uj-limit set of x is closed and 
invariant under the flow of f. If there is a T such that x{t) is contained in a 
compact set for t>T, then the uj-limit set of x is connected. 

Proof See e.g. 0. □ 

The following lemma will be a basic tool in the analysis of the asymptotics, we will 
refer to it as the monotonicity principle. 

Lemma 2. Consider Let U be an open subset of R" and M a closed subset 
which is invariant under the flow of the vectorfield f . Assume there is a function 
F € C(J7, R) such that F(x{t)) is strictly monotone for any solution x(t) of as 
long as x{t) € U H M. Then no solution of whose image is contained inU DM 
has an a- or uj-limit point in U . 

Remark. Observe that one can use M — W\ We will mainly choose M to be the 
closed invariant submanifold of defined by (Q). If one Ni is zero and two are 
non-zero we consider the number of variables to be four etc. 

Proof. Let a; be a solution to (||) which is contained in UDM, and which has maximal 
existence interval (t_, Then F o x is strictly monotone. Suppose p G J7 is an 
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w-limit point of so that there is a sequence i„ — > such that x(tn) — > p. Thus 
F{x{tn)) F{p), but F o x is monotone so that F{x{t)) P{p) sls t ^ — - 
Thus F{q) — F{p) for aU w-hmit points q of x. Since M is closed p € Af. The 
solution a; of (^) with initial value p is contained in Af by the invariance property 
of M, it consists of w-limit points of x so that F{x{t)) — F{p) which is constant. 
Furthermore, on an open set containing zero it takes values in U contradicting the 
assumptions of the lemma. The argument for the a-limit set is similar. □ 

3. Relation between the variables of Wainwright and Hsu and the 

space time metric 

We need to relate the evolution of the variables in the formulation due to Wain- 
wright and Hsu to the evolution of the corresponding spacetime. In order to do this 
we first derive the formulation first presented by Ellis and MacCallum Then 
we derive equations as in [||. Secondly, we turn the argument around and 

start with solutions to the equations of Ellis and MacCallum and then construct 
the spacetime metric. The second part uses the first part. We can then express the 
spacetime metric in terms of the variables of Wainwright and Hsu. Finally, we find 
expressions for the reduced Hamiltonian in terms of the Wainwright Hsu variables. 

3.1. Derivation of the equations of Wainwright and Hsu. Consider first a 
special class of spatially homogeneous four dimensional Lorentz manifolds of the 
form 

(9) {IxG,-dt^ + x^Jit)C^e) 

where / is an open interval, G is a Lie group of class A, Xij is a smooth positive 
definite matrix and the are the duals of a left invariant basis on G. 

We express the condition that a manifold of the form (|^) satisfies Einstein's vacuum 
equations in terms of the variables of Ellis and MacCallum. As these variables are 
defined in terms of a suitable orthonormal basis, we begin by constructing it. Let 
eo = dt and = a^'' Zj, i=l,2,3, be an orthonormal basis, where a is a C°° matrix 
valued function of t and the Zi are the duals of Below, Latin indices will be 
raised and lowered by Sij . 

By the following argument we can assume < Ve(,ei,ej 0. Let the matrix A 
satisfy eo{A) + AB = 0, ^(0) = Id where Bij =< VeoCi, ej > and Id is the 3x3 
identity matrix. Then A is smooth and iS'0(3) valued and if = ^/e^- then 
< Veoe^,e^ >= 0. 

Let 

(10) e{X,Y)=<Vxeo,Y >, 

Oaf} — 0{ea,ef3) and [e^jg^] = 7^560 where Greek indices run from to 3. The 
objects 9af3 and 7^^ will be viewed as smooth functions from / to some suitable M'^ 
and the variables will be defined in terms of them. 

Observe that [Zi, eo] — 0. The span the tangent space of G and < [cq, ej], eo >= 
0. We get 9oo = 9oi ~ and Oa/s symmetric. We also have 7?- = 7oi = and 
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7oj- = We let n be defined as in (||) and 

where we by abuse of notation have written tr(0) as 6. 

We compute the Einstein tensor in terms of n, a and 6. The Jacobi identities for 
ea yield 

(11) eo(ny ) - '^n^^^a^'^ + ^dn,^ = 0. 

The Oi-components of the Einstein equations are equivalent to 

(12) cr^nfcj - n^'^akj = 0. 

Letting bij = 2n^''nkj — tr{n)nij and — bij — ^ti(b)dij the trace free part of the 
ij equations are 

(13) eo(o-ij) + OcTij + Sij = 0. 
The fact that Rqo — yields the Raychaudhuri equation 

(14) eo(0) + %0'^=O 

and using this together with the trace of the ij-equations yields a constraint 

(15) a,,a'' + in,,?^^ - hr{nf) - 



Equations (11 )-([l5[) are special cases of equations given in Ellis and MacCalluni 
|l). At t = we may diagonalize n and a simultaneously since they commute ([12[). 
Rotating Bq, by the corresponding element of 5*0(3) yields upon going through the 
definitions that the new n and a are diagonal at i = 0. Collect the off-diagonal 
terms of n and a in one vector v. By (|l^) and (|l^) there is a time dependent matrix 
C such that v = Cv so that v{t) = OVt since v{0) — 0. Since the rotation was time 
independent < Ve^ei, >= holds in the new basis. 

Introduce, as in Wainwright and Hsu — (Jij/6, Nij — nij/6, Bij — 

2N^'^Nkj — N\Nij and Sij = Bij — \B\5ij and define a new time coordinate 
T, independent of time orientation, by 

Let S+ = |(S22 + S33) and S_ = V3(S22 - S33)/2. If we let N., be the diago- 
nal elements of Nij, equations (|ll|) and (^3|) turn into (^)-(0) The Raychaudhuri 
equation (|lj) takes the form 

(17) e'^-{i + q)e. 



3.2. Construction of the spacetime. Above we derived Einstein's equations 
assuming the spacetime to be of a special form. Next we construct a spacetime of 
the form (H) beginning with initial data as in Definition |^. The construction allows 
us to express the evolution of the metric in terms of the variables of Wainwright 
and Hsu. One step is the following lemma, which also classifies the class A Lie 
algebras. 
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Lemma 3. Ta6Ze|^ constitutes a classification of the class A Lie algebras. Consider 
an arbitrary basis {ci} of the Lie algebra. Then by applying an orthogonal matrix 
to it, we can construct a basis {e[} such that the corresponding n' defined by 
has diagonal elements of one of the types given in table |. 

Proof. Let Ci be a basis for the Lie algebra and n be defined as in (||). If wc change 
the basis according to e- — {A^^)--'ej then n transforms to 

(18) n' ^ {detA)-^A^nA 

Since n is symmetric we assume from here on that the basis is such that it is 
diagonal. The matrix A = diag(l 1 — 1) changes the sign of n. A suitable 
orthogonal matrix performs even permutations of the diagonal. The number of non- 
zero elements on the diagonal is invariant under transformations ( [l8| ) taking one 
diagonal matrix to another. If yl = i^^ij) ^iid the diagonal matrix n' is constructed 
as in ( p^ ) we have n'^.^, = (detA)^^ J2i=i ^ik''^ii so that if all the diagonal elements 
of n have the same sign the same is true for n'. The statements of the lemma follow. 
□ 

Consider class A vacuum initial data {G,g, k) with notation as in Definition |^. We 
construct a spacetime of the form (^) whose induced metric and second fundamental 
form on {0} x G is 5 and k. Let e^, i = 1, 2, 3 be a left invariant orthonormal basis. 
We can assume the corresponding n' to be of one of the forms given in table |^ by 
Lemma ||. The content of (^) is that fc^ = fc(e^,e^) and n' are to commute. We 
may thus also assume kij to be diagonal without changing the earlier conditions 
of the construction. If we let n(0) — n' , 9{0) = — tr^fc and (ly (0) — —hj + 06ij/3 
then (^) is the same as ([1^). Let n, a and 6 satisfy (pi]), (|l3|) and ( p^ ) with initial 
values as specified above. Since (|l^) is satisfied at i = it is satisfied for all times 



For reasons given in connection with ( [15D n and a will remain diagonal so that (12) 
will always hold. 



Let M — I X G where / is the maximal existence interval for solutions to (11)-([L5[) 
with initial data as above. We construct a basis with the same properties as the 
basis used in deriving (pT|)-(p^. Then we define a metric by demanding that the 
basis be orthonormal and show that the corresponding n, tr and 9 coincide with n, 
a and 9. We will thereby have constructed a Lorentz manifold satisfying Einstein's 
vacuum equations with the correct initial conditions. 

Let Hi and ai denote the diagonal elements of n and cr respectively. Let /i(0) = 1 
and fi/ fi = 2ai — 9/3. Let Oi = (Hj-^ifj)^^^ and define Ci = a^e^. Then h associated 
to Ci equals n. We complete the basis by letting eo — df. Define a metric < •, • > 
on M by demanding Cq to be orthonormal with < Bq, Cq —1 and < e^, e,j 1, 
i ~ 1,2,3 and let V be the associated Levi-Civita connection. 

Compute < VeoeijCj >= 0. If 9{X,Y) =< \/xeo,Y > and = 0(e^,ei,), then 
^00 = ^io = 9Qi = 0. Furthermore, 

— eo(aj)(5ij — —9ij 

Oj 

(no summation over j) so that 9ij is diagonal and tr9 = 9. Finally, 
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The constructed Lorentz manifold thus satisfies Einstein's vacuum equations. By 
the above we have. 

Lemma 4. Given initial data as in Definition [J, we get a solution to Einstein's 
equations 

3 

(19) g^^dt' + Y^a-\t)e<E>e 

i=l 

where are the duals of a basis of the Lie algebra of G and 

(20) a,{t) = exp(- / (a, + -e)ds) 



where 9, ai (and Ui) constitute solutions to (|j4)-(|jqj with suitable initial condi- 
tions. If I is the corresponding existence interval, the metric ( [7^ is defined on the 
manifold I x G. 

By Q this metric is gfobally hyperbohc, future causally geodesically complete and 
T — > cx) corresponds to the geodesically complete direction, assuming the initial 
data are not of type IX. 



3.3. The rescaled metric. The equations (19) and ( pOj ) express the spacetime 
metric in terms of the Ellis-MacCallum variables. Let g{t) denote the Riemannian 
metric induced on {t} x G, which we will also consider to be a metric on G. Define 
a rescaled version of this metric by 

(21) m = {a,a^a,f'^{t) ® C- 
We can write it as 

3 

(22) m = Y.^-^m®€ 



where 



Ai(t) — exp( / 2aids). 
Jo 



Defining a new time coordinate by 

(23) 1 = 1 = 

we get 



Jo Jo ^ Jo 



Thus if we consider the to be functions of r we have 

(24) A.(t) =exp( / eE^cir'). 

Jo 

Similarly, 

(25) a,(T)=exp(-/ (3E, + l)ds). 
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However, the can be expressed in terms of S-i- and E_. We have 

(26) El = -h,, . + 1^3 = il^. - 

Thus, if we know what S+ and E_ converge to as t — > cxi, we roughly know how g 
behaves asymptotically. 

3.4. The reduced Hamiltonian. As mentioned in the introduction, there is in 
some situations a subgroup F of the group of diffeomorphisms of G, acting properly 
discontinuously on G, such that g(t) is invariant under T for all t. In such situations, 
we can consider 5 to be a metric on M = G/T. Observe however that the basis e[ 
need not necessarily descend to a basis on M. Let us assume M to be compact. 
We can then consider the reduced Hamiltonian 

Hit) = -[-eit)rv{t) 

where V{t) = vol{M,g{t)). 

Lemma 5. In terms of the time coordinate t defined by we have 

(27) iJ(T) = 6i3(0)F(0)exp(-3 / qds). 



Proof. Since M is compact we can choose a finite partition of unity (f>i i — 1, , 
subordinated to coordinate charts. We can then write 

k 



V{t)^y2f 0» Jdctg{t){dj,d„,)dx\.. dx^ 



JM 

We compute 

dV 



so that 



0V 
dt 



Considering if as a function of r, we get (|7|) using (|^. □ 



4. BlANCHI I 



Bianchi I solutions of (H)-(0) are of the form ((74-, cr_, 0, 0, 0) where (fT+,(T_) are 



independent of time and satisfy 



5. Bianchi II 



Proposition 1. A Bianchi II solution of (Qj-fQj with Ni > and N2 — N3 — 

satisfies 

(28) hm Ni^O 

T — >-00 

and 

(29) lim(S]+,S_) = (a+,a_) 

T — >-00 

where (f7+,f7_) belongs to JCi. 



THE FUTURE ASYMPTOTICS OF BIANCHI VIII VACUUM SOLUTIONS 



11 



Proof. Using the constraint (0) we deduce that 

We wish to apply the monotonicity principle. There are three variables. Let U be 
defined by Ni > 0, M be defined by @) and F(E+,S_,iVi) = E+. Equation ^ 



follows, since if iVi does not converge to zero, we can construct an w-limit point 
in U due to the fact that the variables are contained in a compact set. Combining 
this with the constraint we deduce 

lim q = 2. 

T — >CXD 

Since is monotone and bounded, it converges. Since g — > 2 and the w-limit set 
is connected, c.f. Lemma ^ ( |29| ) holds, but we do not yet know anything about the 
limit. Compute 

/ T. \' 

= 0. 



We get 

I — Zj+ i — (7-1- 

for arbitrary (S-(-,S_) belonging to the solution. Since = (q — 4S+)A^i and 
iVi ^ we have to have cr+ > . If it+ = ^ then cr_ = ± . The two corresponding 
lines in the S]-|-S_-plane obtained by substituting (cr-|-, cr_) into ( |30| ) do not intersect 
any points interior to the Kasner circle. Therefore (T-|- = i is not an allowed limit 
point and the proposition follows. □ 

6. BlANCHI VIq 

Proposition 2. Consider a Bianchi VIq solution of (^-([^ with Ni = 0, N2 > 
and N3 < 0. Then 

(31) lim N2 = 0, lim A^3 = 

r — ^00 r — ^00 

and 

(32) lim E_ = 0, lim Y,+ = -1. 

r— >00 T — *oo 

Proof. By the constraint ^ the variables belong to a compact set. Since 

S'+--(2-g)(l + I]+) 

S+ is strictly monotone as long as q < 2 which is true for the entire solution. The 
monotonicity principle yields 

(33) lim q 2 

T — *00 

for similar reasons as in the proof of Proposition |^. By the constraint we conclude 
that ( ^l|) holds. Since is monotone, ( ^3|) holds and the w-limit set is connected, 
we have 

lim = (cr+,cr-) 
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with a'^ + a'^ — 1. Since N2 and converge to zero and since they satisfy (||) we 
must have 

2{al + (7^) + 2a+ + 2\/3cr_ < 

2(cr2 + ) + 2cr+ - 2\/3o-_ < 0. 

Adding these inequahties we obtain 4 + 4cr+ < 0. Equation ( |3^ ) follows. □ 

Let us now record the fact that Bianchi VIII is very different from Bianchi IX in 
the direction towards the singularity. 

Proposition 3. For every e > 0, there exists Bianchi VIII vacuum initial data 
and a real number > 0, where y^ is closer to the Kasner circle than e, such that 
evolving y^ a time to the past yields a Bianchi VIII point with the property that 

\NiN2\ + \N2N3\ + iNsNil > i. 
Remark. By distance to the Kasner circle we here mean the Euclidean distance in 

Proof. Let e > and x be any Bianchi VIq initial data. By Proposition ^, there is a 
Te > such that applying the flow to x, $(Te, x) is closer to the Kasner circle than 
e/2, where <I>(t, x) represents the solution to with initial data x evaluated in 

T. However, there are Bianchi VIII points as close to ^{T^, x) as we wish. Given an 
77 > 0, we can thus, by the continuity of the flow, choose Bianchi VIII initial data 
with the property that ye is closer to the Kasner circle than e and ^{—Te,ye) is 
closer to x than 77. Choosing x to have E+ = S_ = and N2 — —N3 = and 
letting rj be small enough, we get the conclusion of the proposition. □ 



7. Bianchi VIIq 
The following proposition is a consequence of 

Proposition 4. Consider a Bianchi VIIq solution to ^-Q) with N2 — > 0, 
Ni = and E_ = 0. Then one of the following statements holds 

• = 1 for the entire solution and N2 — — s- cx) as r —> 00. 

• = — 1 and N2 — N3 are constant for the entire solution. 

Given a solution to (||)-(0), if N2{t) ~ N^{t) and S-(t) = for one t then the two 
equalities hold for all r. Thus the above proposition together with the one below 
exhaust the possibilities for Bianchi VIIq solutions. 

Proposition 5. Consider a Bianchi VIIq solution to with N2, N3 > and 

Ni — which never satisfies N2 = and S_ = simultaneously. Then 

lim E+ = -1, lim E_ = 

T — ^00 r — *oo 

and 

lim {N2 - N3) = 0. 

T — ^00 

Remark. In Proposition we also prove that N2 and N3 converge to finite positive 
values. 
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Proof. Assume first that there is an w-Umit point. Then the monotonicity 
principle suffices. Consider 

(34) EV = -(2-g)(S+ + l). 

We apply the monotonicity principle with S_, N3) = S+, M defined by 

(0) and U defined by E?. + {N2 - ^3)^ > and iV2 + A^a > 0. F is strictly monotone 
on a solution in C/n Af since if E'_^{t) = then [iV2 - A^3](r) = so that ^-(t) ^ 
whence [N2 - NsYir) ^ 0, c.f. (|)-(0). If Tfe — > 00 yields an w-limit point, we must 
thus have [S^^ + {N2 - N3)'^]{Tk) ^ or {N2 + N3){Tk) ^ by the monotonicity 
principle. However 

Ij:l + {N2-N^)' 1(1^^1) 

Z_i — — — 



N2N3 N2N3 

has non-positive derivative, and consequently [Et + {N2 — A^3)^](rfc) (if it does 
not converge to zero, Z_i(Tfc) 00). We conclude that E_|_(rfe) —1. By the 
monotonicity of the proposition follows. 

Assume that there is no w-limit point. If there is a sequence ^ 00 such that 
N2{Tk) is bounded, then there is an w-limit point by the constraint Similarly 
if there is a sequence such that N3 is bounded. We can thus assume that and 
A^3 converge to 00. 

The variable S+ is still decreasing and bounded. We can thus assume it to converge 
to (7+ with — 1 < cr_(- < 1. If cr+ = — 1 we are done so assume not. We prove that 
this assumption leads to the consequence that (2 — 17) ^ L^{[0, 00)). Combining this 
with ( |3^ ) we conclude that E+ ^ — 1 contradicting our assumption. 

The intuitive idea is as follows. E+ cr+ but S_ will oscillate between ±(1 — 
roughly speaking. When E_ is small we get a contribution to the integral 
of 2 — g. If E_ spends most of its time close to ±(1 — cr+Y^^ we would have a 
problem, but that turns out not to be the case. 

If N2 and A^3 converge to infinity, we will see that the solution will become oscilla- 
tory. To be more precise, let x and y be defined by 

E_ 



(35) 



(36) y 



(l-S2_)l/2 

V3 N2-N3 



2 (1 - 2^)1/2 • 
Since S^(r) < 1 for all t, x and y are smooth. Let 

(37) 5 = 3(A^2+iV3) + 2(l + S+)xy. 

Then x' — —gy and y' — gx. By the constraint + = 1 so that we can choose 
a 00 such that (x(ro), j/(to)) — (cos(0o), sin(0o))- Define 

(38) = / gis)ds + </)o. 

J To 

Then x{t) = cos(^(t)) and ^(t) = sin(^(r)). 

As mentioned it suffices to prove that 2 — ^^([Oi 00)) under the assumption that 
— 1 < (T-i_ < 1. Consider 

2 - g = 2(1 - - (1 - E^) cos2(^)) = 2(1 - E^) sin2(^) > 2(1 - ct^) sin2(^) 
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where c > by assumption. Consider a period. Assume T is great enough that 
1 < 5(7V2 + A^3)/2 <g< 7{N2 + N3)/2 for t > T and that ^(^2) -<^(ti) = 27r where 
T < Ti < T2. Let us estimate At = T2 — ti . By considering ( p8| ) we conclude that 
At can be chosen arbitrarily small if T is great enough since N2,N3 00. Using 
we derive the consequence that for every e > 

1-e < ; ^ \'; { < l + e 

- {N2 + N3){ti) - 

for every r € [^1,72] if T is great enough. Thus 

At = / dr = / < / — ^ -d^ < 



/?(ri) 5 A(ri) 5(iV2+iV3) " " 5(iV2 + TVs ) (ti ) 
if T is great enough. But then 

/ (2-g)dT> / csin2(^)dT= / csin2(^)-d^ 
/•€(t2) 2 r /■C('^2) 

J5(ri) 7(A/2+A/3) 7(A/2 + A/3)(Ti) 

CTT 5c . , , 

> — At = c'Ar 



7(iV2 + W3)(ri) - 56 

if T is great enough. Here c' > 0. Let M > be any number. Consider an interval 
[to, Tb] where Ta is great enough that the above conditions are met, (Th — Ta)/c' > M 
and ^(rh) — ^(tq) is an integer multiple of 27r. Then 



J (2 - q)dT > M 



and consequently 2 — q ^ L^ii'^j °°))- The proposition follows. □ 
It will be of interest to know that N2 is bounded away from zero. 

Lemma 6. Consider a Bianchi VIIq solution to (^-([^- Then 

N2N3 > c> 

on [0,00), where c is a positive constant. 



Proof. Assume there is a time sequence —^ 00 such that {N2N3)(Tk) converges to 
zero. Since N2 — N3 ^ 0, we conclude that the solution evaluated at converges 
to the Kasner circle. Since the flow is continuous, and Bianchi VIIq belongs to the 
boundary of Bianchi IX, we can construct a sequence of Bianchi IX initial data Xk 
converging to the Kasner circle, such that ^{—Tk,Xk) converges to a Bianchi VIIq 
point with iV2iV3 > 0. This is impossible by Theorem 15.1 of 0. □ 

As was noted in the proof of Proposition |[ the objects x and y defined in (|3^) and 
( ^6| ) can be written as 

x{t) = cos^(t), y(r) = sin^(r) 
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where ^ is defined in ( |38| ) . Since A'2 + N3 is bounded away from zero to the future 
by Lemma ^ and 1 + !]_)_ converges to zero, we can assume that 

(39) g(r) > 2{N2 + N3){t) > c> Vt > tq 

for some positive constant c, where g is defined by (^7|). In other words, the expres- 
sions X and y wiU carry out an infinite number of oscillations as time progresses. 
It is natural to average over a period of this oscillation in order to analyze the 
asymptotics, and we need to know how certain expressions vary over a period to 
be able to do this. 

In the rest of this section we will implicitly assume that all Bianchi VIIq solutions 
satisfy + {N2 - Nsf > 0. 

Lemma 7. Consider a Bianchi VIIq solution to with N2,N3 > 0, and let 

To be such that is fulfilled. Let tq < Ta < Th, and assume 

e(rfc)-e(r,)=7r. 

Then there is a T such that 

(40) ii-i^+^;;-^;i<c, (i+^-^k-^) 



(l + E+)(r2)'- (iV2 + A^3)(w) 
and 

(4,) |l-(^!i±Mll)|<C- + 



(A^2 + ^3)(t2) ' - {N2 + A^3)(W) 

if Ta > T, where Tmax yields the maximum value of N2 + in [Ta,Tb\ and ti,T2 
are arbitrary members of [Ta,rb]. The constant C only depends on the constant c 
appearing in (^3). 



g{s)ds > Cin - Ta), 



Proof. Observe that 



so that 



(42) n-Ta<-. 

c 

Consider 

{N2 + NsY = {q + 2S+)(iV2 + N3) + 2V3E-(^2 - ^3). 

Estimate 

\q + 2E+I < 2|I]+(1 + S+)| + 2J:1 < 2(1 + E+) + 2(1 - E^) < 6(1 + 

and 

|2V3E_(7V2 - iV3)| < 2[S^ + ^{N2- Nsf] < 2(1 - S^) < 4(1 + 
Since N2 + N3 is bounded below by a positive constant for t > tq, we conclude that 

|(iV2 + iV3)'| < Ci(l + E+)(iV2 + iVg) 

for some positive constant C'l . Letting ri , T9 G [tq , Tb] , we get 
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where we have introduced the quantities a and f3 in order to make the foUowing 
arguments easier to fohow. If a > 1, then 



(44) 



< a - 1 < e 



1 



^ kl 

k=l 



If a < 1, then ( |44| ) also holds if we replace a with l/a, since ri and T2 in 
arbitrary. Multiplying this inequality with a then yields the conclusion 

|l-a| <l3e^ 

" and < 1 + 



without any conditions on a. Since we have the bound (n 
we conclude that is bounded by a constant, so that 



are 



< 2, 



(45) 



II- 



(7V2 + iV3)(ri), 



< C2 / (1 + S+)ds. 



(iV2+iV3)(r2) 

Observe that this inequality proves that the left hand side can be chosen to be 
arbitrarily small by demanding to be big enough. Let us estimate the integral 
on the right hand side. Assuming to be great enough that the right hand side of 
( ^5| ) is less than i, and using (^9|) and the fact that S+ is monotone, we conclude 
that 

-d^ < (1 + S4 



(l + E+)ds < (l + E+)(Ta) 



9 



-)ira) 



1 



< 



StT (l + E+)(Ta) 



4 (A^2 + iV3)(w) 

where t^^x corresponds to the maximum value of N2 
this with (^sj) we get 

iN2 



(46) |1 

Consider now 
(47) 
Since 

we conclude that 
(48) |1- 



{N2- 
(1 



(1 



f A^3 in [Ti,r2]. Combining 
)ira) 



(1 



FiV3)(r2)'- "(iV2 + 7V3)( 

f E+)' = -(2-9)(l + S+) 
0<2-g<4(l + S+), 
■5]+)(ri),^^ (l + S+)(ra) 



(1 



E+)(r2) {N2 + N3){ 

by an argument similar to that proving (^6|). By assuming to be big enough 
( ^8| ) leads to the conclusion that we, at the cost of increasing the constants, can 
replace (l + E+)(ro) in the right hand side of ( ^6|) and ( ^8|) with the same expression 
evaluated in t^, yielding the statement of the lemma. □ 

The following lemma constitutes the main observation. 

Lemma 8. Consider a Bianchi VIIq solution to (|^-(0j with N2,N3 > 0, and let 
tq be such that ^3^) is fulfilled. Let tq < Ta < ti,, and assume 
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Then there is aT such that 

(49) I ^ '[E^ - (1 + S+)]ds| <C + ^+fds 

if Ta > T , for some constant C depending only on the constant c appearing in 

Proof. Observe first that 

= (1 - Y.\)x^ = (1 - + = 2(1 + - (1 + Y.+fx^. 

In consequence, the integrand of interest is 

2(1 + S+)a;2 _ (1 + s+) = (1 + S+) cos?7 
where ?7(t) = 2^(r). Let rja = ri{Ta) and rjb — ri{Tb) ^ rja + 2ti and observe that 



(1 + E+) cosTydr = y — cos rjdrj 

l + S+(r;) l + E+(r; + ^ 



, , , cos ridri 
25(77) 25(?7 + ^) ' 

where we in order to obtain the last equality divided the integral into one part 
from rja to 77^ + tt and one from ?7a + to rya + 27r, and then made the substitution 
?7 ^ ry — TT in the second integral. We now wish to prove that the absolute value of 
the integrand in the integral from rja to rja + tt can be bounded by 

[i + Y+in)? . ^ [i + Y+iv + nr 

' 25(77) ^ ' 25(77 + 7r) ' 
since this would prove the assertion. Let us first rewrite 

l + E+(77) l + S+(77 + ^) _ l + S+(77) 5(7;) ^ S+(77 + 7r)-E+(r,) 



25(77) 25(77 + 7r) 25(77) V 9iv + T^)J 25(77 + ^) 

Using ( ^0| ) and the fact that N2 +N3 is bounded from below by a positive constant, 
we note that the second term can be estimated in the desired way; multiply (40) 
by 1 + E+(t2) in order to get an estimate of I]-|_(ti) — Y.+ {t2). If we had 

^ (A^2 + A^3)(^) 

{N2 + N3){rj + T:) 

instead of 

1 _ 

giv + tt) 

in the first term, then we could estimate the first term in the desired way as well. 
However, using (|4^), (^l]) and the fact that 5 and N2 + can be bounded from 
below by a positive constant, one can show that 

(A^2 + ^3)(r?) gjy) 

{N2 + N3){r] + tt) 5(77 + 7r) 
is small enough to yield the desired bound. □ 

Lemma 9. Consider a Bianchi VIIq solution to (j^-(0j with N2,N3, > 0, and let 
tq be such that is fulfilled. Then 

r Y.ids 

(50) Km ; - 1 
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and 

r (2 - q)ds 

J^ j; 2(1.. 

Proof. Observe first that since (|4^) is fulfilled and 1 + S+ ^Oasr^oo, 2 — q ^ 
i^([To,oo)). Moreover, 

2-g<2(l-S^)<4(l + I]+), 
so that 1 + S+ ^ L^{[to, oo))- In order to prove (|50|), let us prove that 

lim — ^ ~ 0. 

Let e > and let T be such that Lemma || can be applied and such that 1 + S+(t) < 
e/(2C) for t >T, where C is the constant appearing in (|49|). Given t > T, let ri 
be the greatest number smaller than r such that C(ti) — ^(T) is an integer multiple 
of TT. We have 

I r[s'_-(i+s+)]dsi< r[E2_-(i+i]+)]dsi 

J Tq J Tf) J T 

+ f - (l + E+)|ds. 

We should divide this expression with (1 + I]_|_)(is, and then evaluate the limit 
as T — > oo. After having carried out this division, the first and the third terms will 
yield something that can be assumed to be arbitrarily small by assuming r to be 
big enough. Using (^), the middle term can be bounded by 

I / '[S^ - (l + S+)]ds| <C f '(l + S+)'ds < C[1 + S+(T)] / \l + E+)ds 
Jt Jt Jt 



<-j^ (1 + 



After having carried out the division, the corresponding contribution can conse- 
quently be assumed to be less than or equal to e/2. Equation ( [50| ) follows. In order 
to prove (^, observe that 

{2-q)- 2(1 + S+) = 2(1 + S+) - 2S^ - 2(1 + 

In consequence the proof of ( |5l| ) is similar to the proof of (|50|). □ 

Lemma 10. Consider a Bianchi VIIq solution to (^-(Qj with 7V2, -^3 > 0, and let 
To be such that is fulfilled. Then, if p > 1, 

1 + S+ e LP([to,oo)). 

Proof. By ( |5T| ) we conclude that for each ai < 2 < 02, there is a T such that t >T 
implies 

ai < ■ 

Introduce the notation 



s 



"^</:(i + E,)ci.<-- 



a = 1 + S+(to) 
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and 



h{T) = I [l + E+(s)]ds. 

J Tn 



By (^) we conclude that 

(52) aexp[—a2h{T)] < h'{T) < aexp[—aih{T)] 

for all T > T. Integrating the leftmost inequality, we get 

— ln[exp[a2^(T)] + aa2{T ~ T)] < h{T) 

for all T >T. Inserting this in the rightmost inequality of (^2|), we conclude that 
h'{T) < a[exp[a2h{T)] + aa2(r - T)]"^ . 

Since h' = 1 + E_|_, and since the quotient ai/a2 can be chosen to be arbitrarily 
close to one, the lemma follows. □ 

Proposition 6. Consider a Bianchi VIIq solution to (^-(j^ with N2,N3 > 0. 
Then 

lim N2(t) = lim N^^t) = uq 

r — ►oo r — ^oo 

where < ng < oo. 

Proof. It is enough to prove that N2N3 converges to a positive number. The reason 
is that this would prove that all the variables are contained in a compact set to the 
future, so that 



N2^{N2N,r/^+N',/' f:r^^ 



1/2 1/2 

converges, since N2N3 converges, TVj + ^3 is bounded from below by a positive 
constant, N2^^ is bounded from above, and N2 ~ hy Proposition ^. Since 

N2- Ns^O, the hmit for TVs follows. Since 

{N2N3y = (2g + 4I]+)iV2TV3, 

the essential point is to prove that 

r 

{q/2 + E+)ds 



To 



converges as r ^ cx). Observe now that 

q/2 + E+ = + E+) + = (1 + - (1 + S+) + 

Since 1 + S+ G L^([ro,oo)) by Lemma |l0|, we need only consider 

^(r) - / - (1 + S]+)]ds. 

J To 

Let T be big enough that (49) is applicable. If T2 > ti > T and T3 is the largest 
number smaller than T2 such that 

^(r3)-^(ri) 

is an integer multiple of vr, then 

Ht2) - V'(ti) = / - (1 + + / - (1 + E+)]ds. 
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The second integral converges to zero as ti and T2 go to infinity, since the integrand 
goes to zero and the length of the interval is bounded (^2|). By ( ^ ) we conclude 
that 

I f + <C f \l + E+)2ds. 

Since 1 + S+ G LP'{[tq, 00)), this expression converges to zero as ti goes to infinity. 
The lemma follows. □ 



8. Bianchi VIII 

We assume that N2, > {) and iVi < 0. Let us begin by giving the asymptotic 
behaviour of the Ni. 

Lemma 11. Consider a Bianchi VIII solution to ^-^) with N2, N3 > and 
Ni <0. Then, 

lim iVi(r) = 0, lim iV2(r) = lim Nsir) = 00. 

r — *oo r — >oo r — >oo 

Proof. Reformulate the constraint as 

+ + ^[N^ + {N2 - mf - 2Ni{N2 + N3)] = 1. 

Observe that all terms appearing on the left hand side are non-negative so that 
the absolute values of, for instance, Ni and N2 — are bounded by numerical 
constants. Assume there is a sequence — > cx) such that A^2('''fe) < C < 00. Then 
the Ni{Tk) are uniformly bounded. By the constraint and E_ are contained 
in a compact set. Thus we can choose a convergent subsequence so that we have 
an a;-limit point ((t+, (t_, ni, ri2, n.3). We apply the monotonicity principle with U 
defined hy N, ^ i ^ 1,2,3, M defined by (|) and TVi, iV2, TVs) = 

N1N2N3. F, evaluated on a solution contained in UnM, is strictly monotone since 
{N1N2N3Y = 3qNiN2N3. If g(r) = then either S'_(t) 7^ or S]'+(t) 7^ by 
(il)"(0)- By the monotonicity principle one of the ni must be zero contradicting the 
growth of \NiN2N3\. We conclude that iV2 ^ 00, but then N3 ^ 00 and TVi ^ 
since Ni{N2 + N3) and A'2 — N3 are bounded by the constraint (^. □ 

8.1. The NUT case. Consider the special case N2 = A'3 and S_ = 0. 

Proposition 7. A Bianchi VIII solution with Ni < 0, satisfying N2 — N3 > 
and I]_ = has the following asymptotic behaviour 

lim E+(r) = lim Ni{t) = 0, lim N2{t) = 00 and lim (iViiV2)(r) = -i 

Proof. The constraint implies 

EV = (l-E^)(l-2E+) + ^iVi^ 

Observe that we cannot have S+ 1 since N1N2 diverges to infinity there, c.f. 
(^. Also, ^ — 1 is an impossibility since S'^ is always positive when S]_|_ is 
close to —1. Since A^i ~> by Lemma pT| and — 1 < S+(t) < 1 Vt by the constraint 
S+ i follows. The last limit follows from the first and the constraint. □ 
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8.2. Oscillatory behaviour. The behaviour of a Bianchi VIII solution as r oo 
is in some sense oscillatory. The quantities that oscillate are E_ and N2 ~ N^. In 
order to analyze the asymptotics, we quantify this behaviour. Let 

(53) i = S_, y= ^(iV2~iV3). 
We have 

i' = -3(iV2 + N3)y + e„ = -(2 - + 3Niy. 

Furthermore, 

y' = 3(iV2 + N3)i + ey, Sy = {q + 2S+)y. 

Observe that by the constraint ^ \x\,\y\,\ex\ and \ey\ are bounded by numerical 
constants whereas N2 + 00 hy Lemma ^ Let g — 3{N2 + 

HI -I 

X = {x, yY and e = {e^, ey)* so that x' ~ Ak + e. Let 

(54) C(t) = / 9{s)ds + 00 

J To 

for some 4>{) and xi{t) — cos(^(r)), yi(r) — sin(^(T)). Then, if xi = (xi,?/i)*, 
x']^ = Axi . Define 



a;i 2/1 



Then = and = 0. Let 

(55) r{T) = {S:^(T)+f{T)f'^ 
and 

(56) x(t) = (a.(T),2/(T))* = (r(ro)cos[C(T)],r(To)sinK(r)])* 

where (/)o has been chosen so that x(ro) = x(to). Since [^(x — x)]' = — <I>e and 
$(t) e 5*0(2) Vt we have 



|i(r)-x(r)||<| r 



.(.)|Ms|. 



By the constraint (|^), we have ||e|| < 9, and thus 
(57) ||i(T)-x(T)|| <9|t-to|. 

The next lemma collects some technical estimates. We present them here in order 
not to interrupt the flow of later proofs. 

Lemma 12. Consider a Bianchi VIII solution of (^-(Qj and let C be a constant. 
Then there is a T depending on C and the initial values such that if [ti , T2] is a 
time interval with ti >T and 

C 



\T2-Ti\ < 



(7V2 + ^3)(t3) 
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for some T3 G [Ti,r2], then 

(58) |S+(fi)-S+(<2)| < 

(59) |[iVi(iV2 + iV3)](ti)- [iVi(Ar2 + Af3)](i2)| < 

(60) |l-(^±^l< 



(A^2 + iV3)(i3) 

C2 

(A^2 + iV3)(t3) 
C3 



{N2 + N3){t2)' - {N2 + N3){t3) 

for arbitrary ti,t2,t^ G [''"ii''2] where Ci, C2 and C3 are constants only depending 
on C . 



Proof. By Lemma 11 we can assume T to be great enough that At = T2 — ti < 1. 
The inequahty 



, {N2 + N^y 
' N2 + N. ' 



< 



foUows from (|)-(0). Thus 



, (^2+^3)(tl) 
\N2+N^){t2) 



I < e«(^^-^i) 



so that if (iV2 
(61) 



iV3)(ii 



II 



> {N2 + N:i){t2) 
{N2 + N^){ti' 

' {N2+N3){t2 



8At 



11 < 8e**Ar 



since Ar < 1. Multiplying this inequality by {N2 + A^3)(t2)/(A^2 + A^3)(ti) < 1 we 
conclude that the assumption {N2 + N3){ti) > (-/V2 + ^3)(i2) is not essential. There 
is thus a constant c < 00 such that 



iN2 + Ns)iti) 

{N2 + N3){t2) 



< c 



for any ti,t2 G [''■i,T2]. Consequently 



(62) 



At < 



C 



< 



(A^2 + iV3)(T3) - (iV2 + iV3)(t3)' 

By (|l)-(0) and |[iVi(A^2 + ^3iri are bounded by constants independent of 

Bianchi VIII solution. Equations ( |5^ ) and ( [59| ) follow from (|6^). Equation (60) 
follows from (|l]) and (||). □ 

Let us give the intuitive idea behind the next lemma. How the variables vary 
during a period is not so interesting, what is interesting is to compute the change 
over a period. In that way the continuous time evolution is replaced with a discrete 
evolution. To make the time step independent of approximation we wish to see 
what happens from x = to the next time i = 0, but then we need to know that 
we have such zeros of i. In order to be able to prove this we assume that r (^5|) is 
not too small. The assumption is not too restrictive; our ultimate goal is to prove 
that r ^ and that a general Bianchi VIII solution in some sense converges to a 
NUT solution. If r is big our iteration is well defined and we will use it to prove 
that r decreases and if r is already small we do not need it. 
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Lemma 13. Consider a Bianchi VIII solution of (^-(Qj- There is a T depending 
on the initial values such that if tq > T and r(To) > [{N2 + iV3)(ro)]~"'^/^, then x 
has at least four zeros Ta < < tq < < Td in 



Taj = [ro - To + 



(iV2 + iV3)(ro)' {N2 + N^){Toy 
such that X is non zero in (Ta,Tb), {Tb,Tc) and {Tc,Td). Furthermore, for any two 
consecutive zeros, for example Ta and Tb, we have 

m> l«r.,-£K)-,|<p_^|^ 

where C is a numerical constant. 

Proof. Let T be great enough that if r e [ti, T2] then 

(64) 3{N2 + N3){t) > 5(iV2 + Ns){to)/2 > 5{N2 + iV3)(r)/4. 

This is possible by Lemma [l^ and Lemma ^ Let x be as in ( |56| ) where (f>Q E [0, 27r] 
in ( |5^ ) has been chosen so that x(ro) = x(ro). If r € [ti, T2] then 

\ar) - e(ro)| = I J^^ g{s)ds\ >^iN2+ N,)iTo)\T - to|. 

Since ^ is a monotonically growing function we conclude 

[0,2^] C [00 - ^^>o + ^tt] C K(Ti),e(r2)]. 

We pick out the zeros in the interval [0,27r]. Let ^(tc,) = 0, S,{Tp) = n/A, ^(r^) = 
Stt/A, ^{ts) — 57r/4 f(Te) — In/A and C(t^) = 2n, where and so on belong to 
[rijTa]. In the intervals [tqjT/j], [t^jT^] and [tcT,^] x cannot be zero if T is great 
enough because of (|5^) and the fact that r(To) > [{N2 + iV3)(To)]"^/^. For the 
same reason x' cannot be zero in [173, r^] and [r^jTe] since x' = —gy + ex. Thus 
there is exactly one zero in [rp, r^] and one in [ts, t^] and no zeros in between. The 
remaining zeros are picked out in the same way. By the above it follows that for 
two consecutive zeros, for example Ta and t^, tt/2 < ^(rt) — ^(tq) < Stt/2. Thus 

- e(ra) -n\< || sin(C(r,) - ^(ra) - n)\ < |(| cos(e(rJ)| + | cos(^(r,))|) 

TT 97r^ 
2r(To) [{N2+ N3) [to 

- [{N2 + N3){Ta)]y^ 

using @ and (||). □ 

Lemma 14. Consider a Bianchi VIII solution of (^-(Qj. There is aT such that if 
Ta and Tb are two consecutive zeros of x and r{Ta) > [{N2 + iV3)(Ta)]-l/2^ Ta > T, 
then if z\ and Z2 are S4- evaluated in Ta and Tf, respectively and wi and W2 are 
~'Ni{N2 + N3) evaluated in Ta and Tf, respectively, we have 

(65) 22 - = [-(1 - Zi)(l + Zi)2 + ^Wi{2 - Zi)](Tb - Ta) + El 

and 

(66) W2- Wi ^ {2zl - 2zi + 2 - 'iWi)wi{Tb - Ta) + £2 
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where 

for i — 1,2 and some numerical constant C . 

Remark. We will prove that — Ta is of the order of magnitude 1/{N2 + N^). Con- 



sequently we can ignore the error terms in ( J6q ) and ( 66 ) as long as the polynomial 
expressions appearing in front of — r^, are of the order of magnitude 1. 

Proof. The idea is to consider the derivatives of and —Ni{N2 + N^) and to 
integrate between Tq and Tb. We will see that and — 7Vi(A^2 + -^3) vary slowly 
so that we can replace them with constants, up to an error we estimate. However, 
E_ and N2 — N3 do not vary slowly so that we will have to use (|7|) and to estimate 
/;>in2(e(r))dr. 

All numerical constants below will be denoted C. Observe that 

[N2 + N3)(Ta) 

by the arguments presented in the proof of Lemma |l^, assuming T to be great 
enough. We will assume T is great enough that ^N2+N3){t) > 5{N2+N3){to)/2 > 
{N2 + N3){t) for all T G [Ta, Tb], as is possible by Lemma ^ We begin by analyzing 
the variation of some relevant expressions in [tq , Tb] . By Lemma we conclude 

m p,M-.,,.„,|<__£_ 

for all T e [Ta,Tf,] and similarly for Ni{N2 + Assume N2 + N3 has a max in 
Tmax in [TajT'b] and a min in Tmm in the same interval. Then 

1 1 



(70) 



' (iV2 + iVa) (w) (A^2 + iVs) ( W) ' 

1 M (A^2 + ^3)(W) , ^ C 



{N2 + A^3)(Tmax) (iV2 + A^3)(rmi„) ' " (A^2 + (Ta) 

by ( |60| ) . Since ^ as defined in (^4|) is a diffeomorphism we can consider functions of 
r as functions of ^. We will be interested in the following integral 

/ sin2 e(r)dr = / sin2(C)-de = — - / sin2(e)dC 



1 



-—- / sin2(e)dC + / sin'iO ( - - ^ He 



TT 



where 



l^i|< 



2ff(Ta) 

c 



Si 



[{N2 + m){TaW/^ 

due to (^3|) and ([70|). Furthermore 

(71) Ta-Tb^ dT^ -di 



5(r„) 9 gi^a) 
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by similar arguments, where 82 is of the same order of magnitude as 5i. Conse- 
quently 

(72) J\m^ aT)dT ^ ^{n - Ta) + 63 

where S3 is of the same order of magnitude as 61. We now have all the necessary 
estimates at our disposal. In [tq , Tt,] we have 

= -(2 - - 35+ = -(2 - 2S^)S+ + 2E^E+ - ^(iV2 - Ns f + 3N^ 



-^Ni{N2 + N3) = -(2 - 2z2)zi + 2zir2(T,) cos^{^) ~ 2r^{Ta) sin^{0 + ^wi + £3 
where 



iN2+N3){ra) 

due to (p9[), the analogous estimate for Ni{N2 + ( p7[ ) and the fact that iVi is 
of the order of magnitude (iV2 + iVs)^^ (^. Integrating, using ([z^), we have 

3 

Z2- zi = [-(2 - 2zi)zi + zir^(ra) - r^(ra) + -wi](Tb - t^) + £4 

where 

k4| < 



(7V2+A^3)3/2(Ta)' 

Using the constraint (|^) to express r^{Ta) in zi and Wi (and A^i), we deduce (65) 
and (|67|), i = 1. The argument to obtain ( |66[ ) is similar. □ 

Lemma 15. Consider a Bianchi VIII solution of (^-(Qj- For all e > there is a 
T such that t >T implies S+(t) < 1/2 + e. 



Proof. Let e > 0. If I]+(t) > 1/2 + e/2 for aU r > T, then -A^i(iV2 + N3) 00 by 
(H)"® ^^'^ Lemma |TT|. Since this cannot occur by (|^) we conclude the existence of 
a.t>T such that Y.+ {t) < 1/2 + e/2. If there is a s > i such that S+(s) > 1/2 + e 
there is an interval [Ti,r2] with ti > t, S]_|_(ti) = 1/2 + e/2, E+(r2) ~ 1/2 + e 
and 1/2 + e/2 < S+(r) < 1/2 + e for aU r e [ti,T2]. Using the constraint (0) to 
eliminate Ni{N2 + N3) in the expression for E'^, we conclude that 

- = I I]V(r)dr < - I N!{r)dr. 

Since A'^{ > — 2eA'^i in [ti, r2] we conclude that 

e ^ 97Vf(n) 
2 - 16e ' 

For T great enough this inequality is impossible by Lemma The lemma follows. 

□ 

Lemma 16. Consider a Bianchi VIII solution of (^-([^- Ni{N2 + N^) does not 
converge to zero. 
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Proof. Assume Ni{N2 + N^) 0. We prove that this assumption forces S+ to 
become negative which, in turn, forces —Ni{N2 + N^) to increase, leading to a 
contradiction. We use Lemma |l^ to achieve the decrease in E+. To apply it we 
need r to be big. However, the constraint yields 

(73) ^ 1 - J:1 + ^NiiN2 + Ns) - ^Nl 

Since the last two terms converge to zero due to our assumption and Lemma 111 
only 1 — is of relevance. From above we have control over S]_|_ by Lemma ]T5| 
and if < — ^ we have nothing to prove as it turns out. 

Let us be more precise. Let T be large enough that S+(t) < | for all t > T. This 
is possible by Lemma Assume also T to be great enough that if t > T and 
-| < 5]+{t), then 

(74) r^t)>l> 



4 - (A^2 + iV3)(i)' 

In order to achieve this we use iM), Lemma pll an d Ni{N2 + N3) ^ 0. Finally let 
T be great enough that Lemma ^and Lemma 1^ are applicable to all tq > T such 
thatr(To)> (iV2+iV3)"'/'(To). 

Let t > T. We prove that we can iterate S+ to become less than or equal to 
— ■=; usmg Lemma Pi If T,+ {t) < -i we are done. The first zero after t, say ti, 
which exists by ([74[) and Lemmajl^, can for the same reason be assumed to satisfy 
S+(ii) > -5. Because of (|6|), (jej), ^) and the fact that iVi(A^2 + N3) 0, we 
conclude 

(75) S+(i2)-S+(ti) < -i(t2-ii) 

if T is great enough and t2 is the first zero after ti. If S]_|_(tj,) > — i, the next 
zero ifc+i will satisfy a relation similar to ([75[). Sooner or later we will thus have 
5j+(ifc) < — ^. Since the distance between two zeros becomes arbitrarily small as 
time goes on and the derivative of S+ is uniformly bounded, cannot at late 
enough times become greater than —j after having been smaller than — i due to 
([75|). Consequently there is a T' such that t > T' imphes ^^(r) < —j. We 
conclude that — iVi(A^2 + A^a) cxd by considering the derivative of Ni{N2 + A3) 
and using Lemma gT|. □ 

Lemma 17. Consider a Bianchi VIII solution to For all e > there is a 

T such that t>T implies -Ni{N2 + A^j) > 1/2 - e. 

Proof. Assuming —A^i(A'2 +A^3) < 1/2 — e/2, we prove that the expression grows at 
late enough times. Under certain conditions on E_|_ this will be seen by considering 
the derivative of —Ni{N2 + N^), but under other other circumstances we apply 
Lemma ^ In order for the iteration to work, we will have to have a good upper 
bound on E+, c.f. (|76|). We will also need Lemma |l^ in order to get a starting 
point greater than or equal to some fixed positive number, irrespective of at how 
late a time we start. 

Let > e > and S be such that 

(76) E+(T)<i + e2 
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for all T > S. Define w = —Ni{N2 + N^)- Since w does not converge to zero there 
is an ?7 > such that for all 5" there is a i > S" satisfying w{t) > rj. Let t>She 
such. We prove that there is &t' >t with the property that 

(77) w{t') >\-\- 
Assume 

(78) w{t) <\-\- 

There are two possibilities. Either E+ < — ^, in which case we will be able to see 
that w increases by considering its derivative. If S+ > — ^ we will be able to apply 
Lemma 

1. There is a T' such that if [ti,t2] is an interval in which < —j and ti > T', 
then 

(79) w{t2)-w{ti)>^witi){t2-ti). 
This can be seen by applying Lemma |ll|, (||) and (|^) . 

2. There is a T" such that if ti < t2 are two consecutive zeros of E_ with T" < ti, 
S+(ti) > -i and 1/2 - e/2 > w{ti) > 77/2, then 

(80) w{t2)~w{ti)>ew{ti){t2~ti). 
We get, in ti 

(81) r^ = l-^l- + ^Ni{N2 + N3) >^e- Se^ 

if r" is great enough. If T" is great enough we can apply Lemma ^ Consider (66) 
with Ta replaced with ti and Tf, replaced with t2- Since 2zf — 2zi + 2 > 3/2 Vzi e M 
and 1/2 — e/2 > w(ti) we get 

W{t2) - W{ti) > ^w{ti){t2 - il) + £2- 

Since M;(ii) > 77/2 we can use @ and (P) to prove that @ holds if T" is big 
enough. This is where we need Lemma |16|. 

Assume t to be greater than T' and T" and that w{t) < 1/2 — e/2 for r e [t, s]. We 
can divide [t,s] into t = so < si < S2 < •■• < Sk < Sk+i — s where I]_|_(si) ~ — i, 
i = 1, fc and S+ is either > — i or < —j in each [s^, s^+i]. If 1 < « < fc then 
we can change Si so that it becomes a zero of I]_ without changing the earlier 
conditions of the construction except that the new I]_|_(si) need not necessarily 
equal — i, c.f. (|8l|). This is due to the fact that we can apply Lemma if t is 
great enough, that | | is uniformly bounded and the fact that the distance to the 
next zero goes to zero as t goes to infinity. 

Since w{t) < 1/2 — e/2 we can apply Lemma |T3to get a zero Si/2 with t < Si/2 
assuming t is great enough and S+(t) > — i, c.f. (|8l|). Since Si/2 — t is of the order of 
magnitude 1/(A'^2 + A''3) and w' is uniformly bounded we can assume w{si/2) > 77/2. 
If S+ < — J in [sg, si], we let = t- Observe that the condition w{v) > r//2 will 
be satisfied for a zero v of S_ with s > v > Si/2 by 1 and 2. If E+ > — ^ in [sk, s] 
we have to take into account the possibility that s need not be a zero of Let 
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in that case Sfe+1/2 be the last zero of S]_ before s. If < —j in [sfc,s] we let 
•Sfe+1/2 = s. We have, applying ( ^ ) or ( fz^ ) in the intervals [si,Si+i] 

w{Sk+l/2) - w(si/2) > ^e?7(Sfc+l/2 - S1/2). 

Using the fact that \w'\ is uniformly bounded and the estimates of Si/2 — t and 
s — Sfc+1/2 we conclude that 

w[s) — w[t) > —€ri(s — t) — 



2 " ' (A^2+iV3)(i) (7V2+iV3)(s) 

for some constant C > 0. We conclude the existence of a t' > t as in @. Finally, 
w cannot become smaller than 1/2 — e once it has become greater than 1/2 — e/2 
given t great enough. We only have to apply the above and observe that \w'\ is 
uniformly bounded and the distance between two zeros of S- becomes arbitrarily 
small as t ^ cxD. □ 

Theorem 1. Let (!]+, !]_, A^i, iV2, he a Bianchi VIII solution to with 
Ni <0 and N2, N3 > 0. Then 



Furthermore 



and 



lim A^i = 0, lim iV2 = 00, lim — 00. 



lim S]_|_ = — , lim S_ = 

T — ^00 2 T — >00 



lim A/'i(iV2 + iV3) = -i, lim (A^2 - iVa) = 0. 



Proof. We begin by proving that for every e > there is a T such that t > T 
implies S+(t) > 1/2 — e. We will then be able to draw the conclusions of the 
theorem. Use the constraint (0) to obtain 

(82) = -^{N2 - iV3)2(I]+ + 1) + ^iVi(iV2 + A^3)(2S+ - 1) + ^NU2 - S+). 

Observe that by Lemma [l7| we can assume —Ni{N2 + N3) to be big. By Lemma [ll| 
we can disregard the term involving Ni. If E+ is smaller than, say, 1/2 — e/2, we 
want E+ to increase. Considering ( |8^ ) two things can happen. Either r is small, in 
which case the first term is negligible and then the second term which is positive 
will dominate. Thus S+ increases. If r is big, we have to use Lemma ^ We make 
these observations more precise in the following two statements. 

1. There is a T' such that if r^(r) < e/36 and S_|_(t) < 1/2 — e/2 for t G [t,s] 
where t > T\ then 

(83) S+(,s)-S+(i)>i(.s-t). 

The proof is as follows. If T' is great enough then —[Ni{N2 + N3)]{t) > ^ for 
r > T' so that if S+(t) < 1/2 - e/2 then 

^7Vi(7V2 + iV3)(2S+-l)>| 

in T. Assume T' to be great enough that we also have 

^[iVf(2-E+)](r)<^ 
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for r > T'. If T e [t, s] we thus conclude that 

Kir) > I 

using (^2|). The statement foUows. 

2. There is a T" such that if r2(r) > e/100 and I;+(t) < 1/2 - e/2 for r G [n, T2] 
where ti > T" anc? ti and T2 are consecutive zeros of E_ then 

(84) E+(r2)-I]+(Ti)> ^(T2-ri) 

The proof is as follows. We have 



(A^2 + iV3)(r) 

for T e [''■1,72] if T" is large enough. If T" is large enough, we can apply Lemma 
p!^ . Consider (^). The relevant polynomial expression is 

-(1 - Zi)(l + Zi)2 + ^11,1 (2 - Zi) = -(1 - Zi)(l + Zi)2 + ^(2 - Zi) 

+ - ^)(2 - ^1) = (^1 + 2)(zi " + liwi - i)(2 - zi). 

Since —1 < zi < 1/2 — e/2 in our situation, and since we can assume T" to be great 
enough that wi > 1/2 - £2/36 by Lemma 0, we conclude 

-(1 - Zl)(l + zi)2 + ^wi(2 - Zl) > ^ - ^ - J. 

Thus (|8^) follows, in which we have absorbed the error term in (^) using (j6^) and 

Let S be greater than T' and T". We start by proving the existence of a i > S" such 
that I]+(t) > 1/2 - e/2. Assume S+(t) < 1/2 - e/2 in [S,s]. We can divide the 
interval by 5* = sq < si < ... < Sk < Sn — s where r^ is either > e/100 or < e/36 in 
each [si, Si+i] and r'^{si) — e/50 for all Si i = 1, k. If S is large enough we can 
assume the Si j = 1, A: to be zeros of I]_ without changing the earlier conditions 
of the construction other than r'^(si) = e/50, i — l,...,fc. The reason is that we 
can apply Lemma |l^ to Si, i — 1, fc, if S is large enough, to get a zero within a 
distance of the order of magnitude 1/(A^2 + -^3)- Using the estimate ( p7| ) and the 
fact that r^ = + , can be made to vary an arbitrarily small amount from the 
original Si to the first zero after it by choosing S large enough. If r'^{T) > e/100 in 
[so, Si], we let Si/2 be the first zero after sq. Otherwise, we let Si/2 — s- We define 
Sfc+1/2 analogously. Using (|3|) and ( p^ ) we conclude 

S+(Sfc+l/2) - S+(si/2) > — (Sfe+1/2 - S1/2) 

assuming e < 10/3. If we use our estimates of Si/2 — S and s — Sfc+1/2 we get 



2 



C C 



10^ ' iN2 + N3)iS) iN2 + N3)is) 

If S" and s — S are large enough this is not possible and there must thus be a t such 
that E+(i) > 1/2 - e/2. 
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If S big is enough, cannot become smaller than 1/2 — e once it has been larger 
than 1/2 — e/2 due to (Q), (p3|), the fact that is uniformly bounded and the 
fact that the distance between two zeros of S_ goes to zero as r ^ oo. 

Combining this observation with Lemma we conclude that 



lim S+(t) = i. 



The constraint (R) then yields 



-N,{N, +Ns)<^{l- \{N, N,f \nI) < ^(1 - Y.\). 

Since the right hand side converges to ^ we can use Lemma |l^ to conclude that 

lim iVi(A^2 + A^s) = -]:■ 

T— >00 Z 

By the constraint we can then conclude that 

lim (iVa - iVs) = 

r — >oo 

and 

lim S_ = 0. 

r— >oo 

□ 



9. Conclusions 

For all class A spacetimes except IX, q converges to a non-zero value as r ^ oo. 
Thus, by (|2^), the following theorem follows. 

Theorem 2. For all Bianchi class A vacuum spacetimes except those of type IX, 
the asymptotic behaviour of the reduced Hamiltonian in the expanding direction is 
given by 

lim i/rcducod(T) = 0. 

r — >oc 

Consider the metric g. By (|2^) we have 

3 

where 

A,(r) = exp( / 6E,dr') 
Jo 

by (p4[). Note that Ej is a vector determined by E+ and 5]_ and that the sum of 
its components is zero. Let us introduce some terminology. 

Definition 4. If two converge to zero and one to infinity, we say that the evo- 
lution exhibits a cigar degeneracy. If one Xi converges to zero and two to infinity, 
we call it a pancake degeneracy. 
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9.1. Bianchi I. In this case the Si are constants satisfying 

3 



2 



^ ' 3' 

i=l 

Thus Xi converges to 0, oo or 1 depending on the solution. We get both pancakes 
and cigars. Concerning the in (25), we observe that if the Bianchi I solution 
corresponds to a special point on the Kasner circle, then two are constant and 
one goes to zero. For all other Bianchi I solutions, one Ui goes to infinity and the 
other two to zero. 



9.2. Bianchi II. The converge to Si where the Si satisfy 

3 „ 

i=l 

The points (si, 32,53) = (f'^i'""!) ones obtained by permutating its 

coordinates are however not allowed as limit points (they correspond to the flat 
Kasner vacuum solutions). Again we get pancakes or cigars depending on the 
solution. As far as the af.s are concerned, one goes to infinity and two go to zero. 

9.3. Bianchi VIo. In this case the Si converge to (|, — i, — i) or one of the points 
obtained by permuting its coordinates. We get a cigar degeneracy. 

Considering the Oi, ( [25| ) shows that one of the Ui goes to zero, but what happens 
to the other two is unclear. 

Proposition 8. Consider a Bianchi VIq solution to Then all the Oi con- 

verge to zero as T —> 00. 



Proof. Considering (^5|) and (|2^), we are interested in the integrals 

/ (l + S+±\/3S_)ds= i / (2 + 2S+±2\/3S_)ds. 
Jo 2 Jq 

Observe however that 

(l + S+)' = -(2-g)(l + E+) 

whence 



l + E+(r) =exph / 
Jo 



(2-<z)ds][l + S+(0)], 



so that, using the expressions for N2 and in (^, we get 

=C2- 



+ 
and 

-N3 

On = 

^ 1 + S+ 

where C2 and C3 are positive constants. Consider now the function 

_ |S2_ + (jV2 - N,)^ 
^ -N2N3 
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This function is monotone decreasing to the future, and bounded from below by 2, 
since 

Thus it converges to a positive real number, and in consequence, the same is true 
of \/Zi. Thus, using the constraint, 

where a is a positive real number. Since N2 and converge to zero as t — > 00 by 
Proposition 0, we conclude that a2 and as converge to zero. □ 



9.4. Bianchi VIIq. There are two possibilities. Either the Si are constants equal 
to (— |, i, i), to (|, — i, — i) or one of the vectors obtained by permuting the coor- 
dinates or they converge to (or a permuted vector). The first possibility 
yields a pancake or cigar degeneracy. The second possibility yields a cigar degen- 
eracy. 

Considering the we have the following proposition. 

Proposition 9. Consider a Bianchi VIIq solution to (^-(0). Then all the Oi con- 
verge to zero as T —> 00. 

Proof. The argument is similar to the Bianchi VIq case, but easier. We just need 
to observe that are bounded away from zero, and 1 + ^ by Lemma 

^ and Proposition ^. □ 

9.5. Bianchi VIII. In this case the converge to (— |, ^, \) up to permutations. 
We have a pancake singularity. 

Let us also consider what happens to the a^ appearing in (|ig|). Considering (p5|), 
two of the Oi converge to zero, but what happens to the third one is more difhcult 
to say, since the integrand tends to zero. Here we wish to prove that it is bounded 
away from infinity. 

Proposition 10. Consider a Bianchi VIII solution to (^-(Qj- The Oi in ( fJ^ j are 

hounded as r —> 00. 



Proof. As already mentioned two of the a,; converge to zero. Let us consider the 
third one. Since Ei — — 2E+/3, we have 



(85) 



ai(r)=exp(/" (2E+ - l)<is) 
Jo 



and in consequence our goal is to prove that E+ does not get that much bigger than 
i . By Theorem |^, we conclude that q — 4E+ — | , so that there is a Ti such that 

Ni{t) < e-^ 

for all T >Ti. Using the constraint (^ to eliminate Ni{N2 + -/V3), we deduce 
S'+ = -(1 - E^ - E2_)(2E+ -l)-^iN,- N, f + ^Nl 
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Observe that in the first term, 1 — S'^ — S'i |. The term involving is 
undesireable, but it can be handled in the following way. Let 

Then 

/' = -(1 - - S2_)(2S+ - 1) - ^(iV2 - Nsf + ^Nf + iq- 4S+)7Vi. 

If S+(r) > i, and if r is greater than some T2, then we can absorb the term 
involving Ni in the term arising from the derivative of iVi , and conclude that 

(86) fir) < -afir), 

where 1 > a > is a suitable constant. If there is a T such that S_|_(t) < ^ for all 
T > T, we are done, so assume not. Let us divide the problem into two subcases. 
Either there is a T > Ti , r2 such that S+ (r) > i for all r > T, or there is a time 
sequence — > 00 such that S+(Tfe) = ^. It is convenient to assume ri > ^i, T2 in 
this case. In the first case, we have 

/(r)</(T)exp[-a(T-r)] 

by integrating (^), and in particular 

0<I]+(T)-i </(r)expha(T-r)] 

so that the integral appearing in ( |85| ) is finite. Consider now the second subcase. 
Assume S+(t) > 1/2 with r > ri. Let i < r be the first point before r at which 
T,+ {t) = \. Then 

/(t) < /(i) cxp[— a(T — t)\ < exp[— i — a{T ~ t)] < exp[— ar] 

since a < 1 and f{t) = Ni(t) < e^*. Thus the integral appearing in ( |85| ) cannot 
diverge to infinity. □ 
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